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Abstract 

• i-H 

j_] ■ The photon structure function F^(x, Q 2 ) has been measured using data taken 

by the OPAL detector at e + e~ centre-of-mass energies of 91 GeV, 183 GeV and 
189 GeV, in Q 2 ranges of 1.5-30.0 GeV 2 (LEP1), and 7.0-30.0 GeV 2 (LEP2), 
probing lower values of x than ever before. Since previous OPAL analyses, new 
Monte Carlo models and new methods, such as multi-variable unfolding, have 
been introduced, reducing significantly the model dependent systematic errors in 
the measurement. 
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1 Introduction 



The measurement of the hadronic photon structure function F^x^Q 2 ) is a classic 
test of QCD predictions The structure function of the photon differs from that 

of the proton because the photon can couple directly to quark charges, as well as 
fluctuate into a hadronic state. The value of F 2 prot °°(i, Q 2 ) exhibits a clear rise towards 
low values of Bjorken x [f|,|5|]. Some theoretical models predict a similar rise in the 
photon structure function |§|-§J, while other models do not require such a rise 0. 
Experimentally, a rise at low-x in the photon structure function has neither been 
observed nor excluded P"U|-P"§|j. The photon structure function is studied at LEP using 
samples of events of the type e + e~ — > e + e~7*7 — > e + e~ + hadrons. The analysis 
presented here uses single-tagged events (from here on referred to as 7*7 events), which 
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means that only one of the scattered beam electrons^ is observed in the detector. These 
events can be regarded as deep inelastic scattering of an electron off a quasi-real target 
photon, and the flux of quasi-real photons can be calculated using the equivalent photon 
approximation [|lj. Figure [I] shows a diagram of this reaction, in which k is the four- 
vector of the incoming electron which radiates the virtual photon, and q and p are the 
four-vectors of the virtual photon and the quasi-real photon, respectively. The symbol 
fq i7 (x,Q 2 ) represents the parton densities of the quasi-real photon. 




e anfitag (P ~ 0) 



Figure 1: Deep inelastic electron-photon scattering. 



The cross-section for deep inelastic electron-photon scattering can be written in 
terms of structure functions as IH1 



d 2 Cr e7 ^eX _ 27TQ 2 

dxdQ 2 xQ 4 



(l + (l- y f)F2(x,Q 2 )-y 2 FZ(x,Q< 



where Q 2 = —q 2 is the negative value of the four-momentum squared of the virtual 
probe photon, a is the fine structure constant and x and y are the usual dimensionless 
deep inelastic scattering variables, defined by 

I = ^1, B ™. (2) 
2p ■ q p ■ k 

In the kinematic region studied in this paper, )/ 2 < 1, so the contribution from the 
longitudinal photon structure function in Equation [l] is neglected. 

In contrast to measurements of the proton structure function, here the energy of the 
target particle is not known. In consequence, the kinematics cannot be fully determined 
without measuring the hadronic final state, which is only partially observed in the 

1 For conciseness, positrons are also referred to as electrons. 
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detector. This leads to a dependence of the measurement on Monte Carlo modelling 
of the hadronic final state, which enters when an unfolding process is used to relate 
the visible distributions to the underlying x distribution. 

The analysis presented here uses OPAL data collected during the years 1993-1995, 
1997 and 1998, at e+e" centre-of-mass energies of 91 GeV (LEP1), and 183 GeV 
and 189 GeV (LEP2). This is the first OPAL measurement of the four-flavour photon 
structure function using the 183 GeV and 189 GeV data, though the charm contribution 
to F2 has already been measured [GO]. In this analysis new Monte Carlo programs 
and improved unfolding methods have been introduced, which are also used in the 
re-analysis of the LEP1 data. 



2 The OPAL detector 



The OPAL detector is described in detail elsewhere [pi]] ; only the subdetectors which 



are most relevant for this analysis, namely the low angle electromagnetic calorimeters 
and the tracking devices, are detailed below. The OPAL detector has a uniform mag- 
netic field of 0.435 T along the beam direction throughout the central tracking region, 
with electromagnetic and hadronic calorimetry and muon chambers outside the coil. 

The forward detectors (FD) cover the 9 region from 60 to 140 mrad at each end of 
the OPAL detector^. They consist of cylindrical lead-scintillator calorimeters with a 
depth of 24 radiation lengths (Xq) divided azimuthally into 16 segments. The energy 
resolution for electromagnetic showers is 18%/ yE, where E is in GeV. An array 
of three planes of proportional tubes buried in the calorimeter at a depth of 4 X 
provides a precise shower position measurement, with a typical resolution of 3-4 mm, 
corresponding to 2.5 mrad in 9, and less than 3.5 mrad in <p. 

The small- angle silicon tungsten lumino meter (SW) covered the region in 9 from 
25 to 60 mrad from 1993-1995. For LEP2 running, a radiation shield was installed 
which moved the lower edge of the useful SW acceptance to 33 mrad. The SW detector 
contains 19 layers of silicon alternating with tungsten. Each of the 16 azimuthal wedges 
is divided into 64 pads for positional measurement. The energy resolution of the SW 
detector is about 24%/a/E at LEP1 energies, and about 6% at LEP2. 

Charged particles are detected by a silicon microvertex detector, a drift chamber 
vertex detector, and a jet chamber. Outside the jet chamber, but still inside the 
magnetic field, lies a layer of drift chambers whose purpose is to improve the track 
reconstruction in the ^-coordinate. The resolution of the transverse^] momentum for 



2 In the OPAL right-handed coordinate system the x-axis points towards the centre of the LEP 
ring, the y-axis points upwards and the z-axis points in the direction of the electron beam. The polar 
angle 9 and the azimuthal angle <j> are defined with respect to the z-axis and x-axis, respectively. 

3 Transverse is always defined with respect to the z-axis. 
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charged particles is <J p Jpt = y0.02 2 + (0.0015p t ) 2 for | cos 6 1 ! < 0.7 , where p t is in 
GeV, and degrades for higher values of | cos#|. Tracks are accepted up to a limit of 
| cos#| < 0.9622. 

Both ends of the OPAL detector are equipped with electromagnetic endcap calorime- 
ters covering the range from 200 to 630 mrad in polar angle. They are homogeneous 
devices composed of arrays of lead-glass blocks of 9.2 x 9.2 cm 2 cross-section and typ- 
ically 22 X in depth, giving good shower containment. In the central region, outside 
the solenoid, is the electromagnetic barrel calorimeter of similar construction. 

The deep inelastic scattering events are triggered with high efficiency by the large 
energy deposits of the scattered electron in the forward calorimeters (FD and SW) and 
by charged particles seen in the tracking devices. 



3 Kinematics and data selection 



To measure F^x, Q 2 ), the distribution of events in x and Q 2 is needed. These variables, 
illustrated in Figure |], are related to the experimentally measurable quantities 

Q 2 = 2£ b £ tag (l-cos# tag ) (3) 

and 2 

X = Q 2 + W 2 + P 2 (4) 

where Ei, is the energy of the beam electrons, E tag and 9 ta , g are the energy and polar 
angle of the deeply inelastically scattered (or 'tagged') electron, W 2 is the invariant 
mass squared of the hadronic final state and P 2 = —p 2 is the negative value of the vir- 
tually squared of the quasi-real photon. The requirement that the associated electron 
is not seen in the detector ensures that P 2 < Q 2 , so P 2 is neglected when calculating 
x from Equation f|. The electron mass is neglected throughout. 

Three samples of events are studied in this analysis, classified according to the 
centre-of-mass energy and the subdetector in which the tagged electron is observed. 
Data from LEP1 are used, with e + e~ centre-of-mass energies close to the Z° mass, as 
well as LEP2 data with centre-of-mass energies of 183 GeV and 189 GeV. Electrons are 
tagged using the SW detector at all centre-of-mass energies, since accessing the lowest 
possible x region requires measuring the electrons scattered at the lowest possible 
angles. The FD is only used for the LEP1 data, to provide a sample in the same range 
of Q 2 as the LEP2 sample, but using a different subdetector for tagging the scattered 
electron. The three samples used are termed LEP1 SW, LEP1 FD and LEP2 SW. 
Each sample is further split into two bins of Q 2 . 

Events are selected by applying cuts on the scattered electrons and on the hadronic 
final state. The cuts are listed in Table [I]. A tagged electron within the clear acceptance 
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of SW or FD is selected by requiring £ tag > 0.75 E h at LEP1 or E tag > 0.775 E h at 
LEP2. This cut effectively eliminates events originating from random coincidences 
between off-moment umQ beam electrons faking a tagged electron and untagged 77 
events. The cut is higher at LEP2 because of the larger off-momentum background 
in the 183 GeV data. To ensure that the virtuality of the target photon is small, the 
highest energy cluster in the hemisphere opposite the tagged electron must have an 
energy E a < 0.25 E\> (the anti-tag condition). To reject background from 7*7 events 
with leptonic final states, the number of tracks in the event passing quality cuts and 
originating from the hadronic final state, iV c h, must be at least three, of which at least 
two tracks must not be identified as electrons. The tracks and the calorimeter clusters 
are reconstructed using standard OPAL techniques |22| which avoid double counting 



of particles which produce both tracks and clusters. Finally, the visible invariant mass 
W v i S of the hadronic system based on tracks and calorimeter clusters and including 
the contribution from energy measured in the forward calorimeters (FD and SW) is 
required to be in the range 2.5 GeV < W w - ls < 40 GeV for the LEP1 sample and 
2.5 GeV < Wvi B < 60 GeV for the LEP2 sample. The distribution of W vis extends to 
higher values at LEP2 because more energy is available from the beam electrons. In 
addition, the background from Z° — > hadrons, which the maximum W Y - 1S cut is designed 
to reject, is lower at LEP2 than at LEP1. 

The trigger efficiencies were evaluated from the data using sets of separate triggers 
and found to be larger than 99% for all of the samples. 

The cuts applied to each sample are listed in Table |I| and shown as dotted lines in 
Figures The numbers of events in each sample passing the cuts, the integrated 
luminosities and the Q 2 ranges are listed in Table |2|. The luminosity for the LEP1 data 
is lower for the SW sample than for the FD sample because the SW sample does not 
include data from 1995. 



cut 


LEP1 SW 


LEP1 FD 


LEP2 SW 


E tSug /E h min 


0.75 


0.775 


#tag m i n [mrad] 


27 


60 


33.25 


# tag max [mrad] 


55 


120 


55 


EJ E h max 


0.25 


N ch min 


3 (2 non-electron tracks) 


W vis min [GeV] 


2.5 


Wyis max [GeV] 


40 


60 



Table 1: The selection cuts applied to each data sample. 



4 Off- momentum electrons originate from beam gas interactions far from the OPAL interaction 
region and are deflected into the detector by the focusing quadrupoles. They appear predominantly 
in the plane of the accelerator. 
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(Q 2 ) [GeV 2 ] 


sample 


luminosity [pb 1 




events 


Q 2 range [GeV 2 ] 


1.9 


LEP1 SW 


74.6 


4356 


1.5-2.5 


3.7 






4010 


2.5-6.0 


8.9 


LEP1 FD 


97.8 


1909 


6.0-12.0 


17.5 






1578 


12.0-30.0 


10.7 


LEP2 SW 


222.9 


4593 


7.0-13.0 


17.8 






5495 


13.0-30.0 



Table 2: The integrated luminosity, number of selected events, and Q 2 range for each 
data sample. The error on the luminosity is negligible. 



4 Monte Carlo and background 



Monte Carlo programs are used to simulate 7*7 events and to provide background esti- 



mates. All Monte Carlo events are passed through the OPAL detector simulation |23 
and the same reconstruction and analysis chain as real events. 



It has been seen in previous studies [^,|T^,|T9[ that Monte Carlo modelling of the 
hadronic final state is a large source of systematic error in the F% measurement. While 
there is still no 7*7 Monte Carlo generator which describes all the features of the 
observed hadronic final state, improved Monte Carlo programs have become available 
since the previous OPAL measurements of F^ and it is now possible to reject some of 
the models which do not describe the data adequately. 

The Monte Carlo generators used to simulate signal events are HERWIG 5.9 |2"4 



HERWIG 5.9+fc t (dyn) J25f, PHOJET 1.05 g§], and F2GEN |7|. HERWIG 5.9 is a 



general purpose Monte Carlo program which includes deep inelastic electron-photon 
scattering. The HERWIG 5.9+A;t(dyn) version uses a modified transverse momentum 
distribution, k t , for the quarks inside the photon, with the upper limit dynamically 
(dyn) adjusted according to the hardest scale in the event, which is of order Q 2 . The 
distribution was originally tuned for photoproduction events at HERA ||28|| . PHO- 
JET 1.05 simulates hard interactions through perturbative QCD and soft interactions 
through Regge phenomenology and thus is expected to provide a more complete pic- 
ture of 7*7 collisions than HERWIG 5.9. It is used here for the first time in an OPAL 
F2 analysis. In F2GEN, 7*7 events are generated with a pure qq final state. No QCD 
effects between the two quarks are simulated, and the pointlike mode was used, in 
which the angular distribution of the qq final state is taken to be the same as for a 
lepton pair. The generated luminosities of the Monte Carlo samples were 3-6 times 
the data luminosity. 

A sample was generated with each of the Monte Carlo models using the GRV 
LO parameterisation of F%, taken from the PDFLIB library ||29|| , as the input 
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structure function. This version assumes massless charm quarks. To study the effect 
of a different input structure function on the final state, another sample was generated 
using HERWIG 5.9Q with the SaSID J7[ parameterisation of F 2 7 . 

All of the Monte Carlo programs except PHOJET 1.05 allow generation of events 
and cross-section calculation according to a chosen structure function. PHOJET 1.05 
uses an input structure function for the total cross-section calculation but always pro- 
duces the same x and Q 2 distributions. The measurement of F% requires knowledge 
of the structure function used to generate these distributions, so the effective inter- 
nal structure function in PHOJET 1.05 was found by comparing the generator level 
distributions to those of HERWIG 5.9+fc t (dyn) with the GRV LO structure function. 
The ratio of the x distributions in the two samples for each Q 2 range was fitted to 
a polynomial, which, after multiplying by the GRV LO structure function, gave the 
structure function for the PHOJET 1.05 sample. 



The dominant background comes from the reaction e + e~ — > e + e~7*7 — > e + e~T + T~ 
proceeding via the multiperipheral process shown in Figure [I], with e + e _ — > e + e~7*7 — > 
e + e~e + e~ giving a smaller contribution. These were simulated using the Vermaseren 
The process Z° — > hadrons also contributes significantly at LEP1 and was 
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program 

simulated using JETSET [[31 1. Because the aim is to measure the structure function 
of the quasi-real photon, the hadronic 7*7* events are also treated as background. 
They were generated using PHOJET 1.05 with the virtuality of the target photon 
restricted to P 2 >1.0 GeV 2 (LEP1) or 4.5 GeV 2 (LEP2). Other sources of background 
considered were Z° — ► t + t~, simulated with KORALZ ||32|| , non- multiperipheral four- 
fermion events with e + e~r + r~ and e + e~qq final states, which were simulated with 
FERMISV 1 53]|, and W + W~ — > hadrons and untagged 77 events, which were both 
simulated with PYTHIA [[34}| . The contributions from all these were found to be 
negligible in all the samples. The number of expected signal and background events 
for each data sample is shown in Table EjL 



(Q 2 ) [GeV 2 ] 


signal 


7*7 — > T + T 


7*7 — > e + e 


Z° -> had. 


7*7* — > had. 


1.9 


4237±66 


85±6 


7±1 


10±2 


18±3 


3.7 


3834±64 


115±7 


9±1 


18±3 


33±4 


8.9 


1749±45 


85±6 


3±1 


42±6 


30±4 


17.5 


1378±41 


71±6 


4±1 


84±8 


42±5 


10.7 


4403±68 


135±6 


4±1 


9±1 


25±2 


17.8 


5171±74 


220±7 


6±1 


14±1 


45±2 



Table 3: The measured number of signal events (observed events corrected for back- 
ground) in each data sample, and the expected number of background events according 
to Monte Carlo. Statistical errors only are given. 



3 The fct(dyn) version was not used for this test. 
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Figures |2|, and ^ show comparisons between data and Monte Carlo distributions 
for variables relating to the scattered electrons. The main differences are in the LEP1 
SW data sample, where the observed cross-section for selected events is significantly 
higher than either of the Monte Carlo predictions. The LEP1 SW data sample is more 
peaked than the HERWIG 5.9+A;t(dyn) prediction at low # tag , and there are large 
discrepancies between the LEP1 SW data sample and the Monte Carlo samples in the 
distribution of E^/E^ (Figure ||d), with the data having an excess at the low end of 
the plot. However, in this region, the anti-tag distribution is more influenced by the 
hadronic final state than by the anti-tagged electrons, which are almost all above the 
cut. 



In the distributions of variables related to the hadronic final state (Figures |5HTTD, 
there are large discrepancies. Figure § shows the hadronic energy flow as a function 
of pseudorapidity. It can be seen that HERWIG 5.9 tends to put too little energy in 
the central region of the detector, while F2GEN has a much higher peak in the central 
region than the data. PHOJET 1.05 gives the best description of the data in these 
plots. 

The large differences between HERWIG 5.9+A; t (dyn) and PHOJET 1.05 in the W vis 
distributions (Figures |a, |]a and |7]a) and the x v ; s distributions (Figures ||b, ||b and ^b) 
are mostly due to the fact that PHOJET 1.05 does not use the input structure function 
to generate the x distribution, although there are also differences arising from the final 
state modelling. The quantity x vis is calculated by inserting W vis into equation (|. For 
the number of tracks, N cil , in the LEP1 SW sample (Figure |]c), the data distribution 
is above the Monte Carlo distributions at high N^. The Monte Carlo descriptions of 
iVch are closer to the data in the LEP1 FD and LEP2 SW samples (Figures ^|c and 0c). 
There are significant differences in the distributions of E™ 1 (Figures [|d, |6]d and [7|d on 
a linear scale and Figures |S], |It] and [TTJ on a logarithmic scale and divided into three 
bins of x v j S ). E^ is the transverse hadronic energy out of the plane containing the 
beam line and the tagged electron: 

£° ut = E^,t|sin^| (5) 

i 

where i runs over all particles in the hadronic final state, with transverse energy E^. 
The azimuthal angle 0, is measured from the tagged electron. HERWIG 5.9 has too 
few events at large E™ 1 and F2GEN has too many, especially at low values of a; v i S . 



Also shown in Figures |9|, [T(] and [TT] is E ior /E tot , the observed energy in the forward 
region divided by the total observed energy. Of the four Monte Carlo models, PHO- 
JET 1.05 shows the best agreement with the data in the hadronic final state, HER- 
WIG 5.9+A; t (dyn) describes the data better than HERWIG 5.9, and F2GEN gives the 
worst description of the data. For the final results, only PHOJET 1.05 and HERWIG 
5.9+/c t (dyn) are used; however, all four Monte Carlo models are used in studies of the 
unfolding methods. 

None of the Monte Carlo programs discussed above contain radiative corrections 
to the deep inelastic scattering process. These are dominated by initial state radiation 
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from the deeply inelastically scattered electron. The final state radiation is experimen- 
tally integrated out due to the finite granularity of the calorimeters. The Compton 
scattering process contributes very little, and the radiative corrections due to radiation 
of photons from the electron that produced the quasi-real target photon were shown 
to be very small [35|. 



In this analysis the radiative corrections have been evaluated using the RADEG 
program |]36f| , which includes initial state radiation and the Compton scattering pro- 
cess. The calculations are performed using mixed variables, which means that Q 2 is 
calculated from electron variables, while W 2 is calculated from hadronic variables. The 
value of x is found from W 2 and Q 2 , exactly as for the experimental analysis. The 
RADEG program allows calculation of differential cross-sections in bins of x and Q 2 , 
with or without radiative corrections, while applying the experimental restrictions on 
the electron energy and polar angle, the minimum invariant mass and the anti-tag 
angle. The predicted ratio of the differential cross-section for each bin of the analysis 
is used to correct the data, i.e. the measured is multiplied by the ratio of the 
non-radiative and the radiative cross-sections. The radiative corrections reduce the 
cross-section in the phase space of the present analysis. They are largest at low values 
of x. The size of the largest correction to any bin is 14%. The radiative corrections are 
insensitive to the choice of F% in the calculation, for example, the difference between 
the predicted corrections for the GRV and the QPM structure functions is negligible 
in all regions of x and Q 2 . 

Several theoretical ansatzes exist for how F% should behave as a function of P 2 |37} - 



39]. They all predict a decrease of F^{x, Q 2 , P 2 7^ 0) with P 2 which is strongest at low 
values of x. This means that applying a correction to obtain an F% which is valid at 
P 2 = would change the shape of the measured F$- The size of the effect was studied 
based on the GRS |J7[] and SaSID |SS| parameterisations using three values of P 2 - 



0.01, 0.02 and 0.03 GeV 2 , and four values of Q 2 : 1.9, 3.7, 9.5 and 17.6 GeV 2 . These val- 
ues reflect the range of (P 2 ) values expected in the data, and the (Q 2 ) values of the data 
samples. The largest suppression observed at x — 0.003 for P 2 = 0.01/0.02/0.03 GeV 2 
is around 7/12/17%; however, at low x the predictions differ by more than a factor 
of two. Because the distribution of P 2 in the data is not known and the theoretical 
predictions differ significantly, no correction is applied. 



5 Determination of _F 2 



5.1 Unfolding 

Unfolding is a statistical technique which is used in this analysis to correct the measured 
distributions for detector effects. The unfolding problem can be described as follows: 
A measurement gj is made of a true distribution /j, with j and i bins respectively. The 



11 



gj distribution differs from /, because of 



• Limited resolution: 

Events in a certain bin of the true distribution can be spread into several bins 
in the measured distribution. In 7*7 events this can be due to missing energy 
in the forward region (FD and SW plus the beampipe), as well as the intrinsic 
energy resolution of the subdetectors. 

• Limited acceptance: 

Some events fail the cuts and therefore contribute only to the true x distribution, 
and not to the measured distribution. 

• Background: 

There are also events originating from other reactions which only contribute to 
the measured distribution. 



Resolution and acceptance effects can be described in terms of a response matrix Ay. 
The distributions gj and /, are then related by 

n 

9 3 Y, A :if> ■ Ij j ( 6 ) 

i=l 

where bj is the background distribution. Equation |6] cannot simply be inverted to find 
fi because this does not lead to a stable solution^. Instead, after bj is subtracted, 
fi is estimated from a procedure which includes a smoothing (or regularisation) of 
the distribution to reduce statistical fluctuations. The response matrix Aij is found 
from Monte Carlo simulation. After unfolding, F% is determined by reweighting the 
input structure function of the Monte Carlo according to the ratio of the unfolded x 
distribution to the x distribution in the Monte Carlo: 

where -F^mc * s the input structure function of the Monte Carlo sample used for unfold- 
ing, A^Mc.i is the number of events in the Monte Carlo sample in the zth bin and f(x) 
is the unfolded x distribution, obtained from a quadratic spline fit to fi, with f\x) = 
at the edges of the acceptance region. represents the measured value of the average 
K7 of the events in bin i of the unfolded distribution, and not the average over x. The 
distinction is important mainly at low x, where f(x) changes most rapidly. 



The main unfolding program used for this analysis was GURU the programs 



RUN JO]] and BAYES [J2| were used to check the unfolding. In previous OPAL F2 
analyses, only RUN was used. The programs GURU and RUN are similar in principle, 
but have different implementations. Due to its use of standard matrix techniques, 



See the GURU or RUN Q] documentation for an explanation of this effect. 
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GURU has a much simpler structure and is easier to modify than RUN, which uses a 
maximum likelihood fit with spline functions. The BAYES program performs unfolding 
using an iterative method based on Bayes' theorem, leaving regularisation to the user 
(e.g. a polynomial fit after each iteration except the last). The programs GURU, RUN 
and BAYES were compared by unfolding test distributions including experimental data. 
The result of one such study, unfolding the LEP1 SW data with HERWIG 5.9+/c t (dyn), 



using the measured distribution x v i S , is shown in Figure |12|a. The results from the three 
programs are generally consistent, though the BAYES program tends to assign smaller 
errors than the other two. 

The amount of regularisation in the unfolding is set by the number of degrees of 



freedom (NDF). It can be seen in Figure |T2| b that increasing this number leads to larger 



unfolding errors. Reducing the number of degrees of freedom increases the correlations 
between the bins. The number of degrees of freedom to be used is determined by the 
GURU program from statistical analysis of the data distributions, as described in the 
GURU documentation |40| . 



The unfolding of the variable x is done on a logarithmic scale in order to study the 
low-x region in more detail. Several improvements to the unfolding procedure used in 
previous OPAL analyses have been implemented. These are described in the following 
sections. 



5.2 Reconstruction of W 



The visible hadronic mass squared is defined by 

W3. = (E^) 2 - (E^) 2 - (E^) 2 - (E^) 2 (8) 

i i i i 

where % runs over all tracks and clusters in the hadronic final state, with energy Ei 
and momentum vector (pi, x ,Pi,y,Pi,z)- Because hadronic showers are not well contained 
in the electromagnetic calorimeters, and some energy is lost in the beampipe, the 
measured energy in the forward region is less than the true energy. It is possible to 
improve the reconstruction of W v - ls by including kinematic information from the tagged 



electron @. W v 2 is can also be written as 

w^ = (£ft+)Eft-)-Eft,t) a ( 9 ) 



where pi± = Ei ± p^ z and pi tt = \l Pi, x 2 + Pi, y 2 ■ Using conservation of energy and 
momentum, and assuming that the untagged electron travels along the beam direction, 
the sums over pi + and p^t can be replaced by electron variables, leading to 

Wrec = (Pbeam+ ~ Ptag+) (E Pi~) ~ (Ptag.t) 2 (10) 

i 

where Pbeam+ and p tag + are calculated for the tagged electron before and after scatter- 
ing, respectively, and pta g ,t is calculated according to the definition of p i>t , above, for 
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the tagged electron. When using Equation |TU| instead of Equation [5] to evaluate the 
measured hadronic mass of each event, the hadronic energy enters only through the 
Pj_ term. This is an advantage because the leptonic energy resolution is usually better 
than the hadronic energy resolution. The (reconstructed) variable formed in this way 
is called W Tec . The corresponding measurement of x is called x rec . The improvement 
in the resolution when using W Tec instead of W V i S depends on how well the hadronic 
system is measured, and is therefore both detector dependent and model dependent. 

Even after the above technique has been applied, the value of W Tec is still gener- 
ally smaller than the true value. This is mainly due to energy losses in the forward 
calorimeters, in which only about 40% of the hadronic energy is observed. This makes 
the measurement very dependent on the angular distribution of the hadrons in the final 
state. In an attempt to make the energy response of the detector more uniform and 
to reduce the systematic error due to the uncertainty in the Monte Carlo modelling, a 
new (corrected) variable is formed: W COT , along with the corresponding x measurement 
x COT , in which the contribution from the forward region has been scaled by a factor of 
2.5. This factor was obtained by comparing the generated and measured energy in the 
forward region in Monte Carlo events. The uncertainty of this factor has been taken 
into account in the evaluation of the systematic errors. 

Wcor = (Pb e am+ ~ Ptag+ ) Pi- ) ~ (Ptag.t ) 2 ( 1 1 ) 



for clusters observed in SW or FD 
l ~ I p,-_ otherwise 




Figure |13| a shows the correlation between the generated W and the three measured 
quantities W v - m , W Tec and W cor for the LEP2 SW sample generated with HERWIG 
5.9+fct(dyn). The spread of W cor in W is larger than that of the other two variables, 
because in W cor more weight is given to the forward region, where the energy resolution 
is worse than in the central region. In general, W COT is still lower than W, mainly 



because of energy lost in the beampipe. Figure |13| b shows the correlation between the 
generated energy in the forward region and the scaled observed energy in that region, 



y cor • 



The three measured variables x v i S , x rec and x COT can be used to unfold the true 
variable x. The difference between the results is model dependent and generally small 
when using Monte Carlo models that already give a good description of the hadronic 
final state (Figure pT2|c) . 



5.3 Two dimensional unfolding 

The GURU program can also be used to perform unfolding in two dimensions. As with 
the attempts to improve the W reconstruction described above, the motivation is to 
reduce the dependence of the unfolding on a particular Monte Carlo model. There is 
information in every event about the angular distribution of energy in the detector, but 
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it is not fully exploited if only x is used in the unfolding procedure. Including another 
variable in a second unfolding dimension makes more direct use of this information. 

Two variables were considered as possible second unfolding variables. They were 
chosen because they are very sensitive to the angular distribution of the hadrons in the 
final state: 



• E^/Etot, the transverse hadronic energy out of the plane containing the beam 
line and the tagged electron (see Equation |), divided by the total observed 
energy. 

• Ef or /E tot , the observed energy in the forward calorimeters divided by the total 
observed energy. 



These variables are shown in Figures |9HTl|. To test the two dimensional unfolding 
procedure, a random number was used as the second variable. The results are shown 
in Figure |T2|d. They are consistent with one dimensional unfolding, which is as expected 
since there is no extra information in a random number. 



5.4 Testing the unfolding methods 

To find which of the unfolding methods was the most reliable, OPAL data was un- 
folded with four Monte Carlo programs HERWIG 5.9, HERWIG 5.9+fc t (dyn) PHO- 
JET 1.05 and F2GEN. Despite not giving good descriptions of the data, HERWIG 
5.9 and F2GEN were included in this study to investigate the effectiveness of the new 
techniques using extreme models. The best methods are considered to be those which 
give the smallest difference between the unfolded results with the four Monte Carlo 
models. The quantity compared is x 2 , defined by 



r = iEE i^P 1 ) (12) 

models * ' ' 



where is the value of the unfolded result in the zth bin, (F^) is the average of 
the results from all four Monte Carlo models and <7j is the statistical error for F^. 
The values of \ 2 are shown in Table ^ and the results of the unfolding are shown in 
Figures P^-fT5|. The \ 2 values show how large the differences between the models are, 
compared to the statistical error. In general, the lowest values of x 2 were obtained 
using two-dimensional unfolding with x cor as the first variable and E™ 1 / E tot as the 
second variable, which consequently is used as the standard unfolding method for the 
results. The \ 2 values are usually smaller for the LEP1 FD sample than the other two. 
This is partly due to the smaller number of bins. Additionally, the lower statistics 
in the LEP1 FD sample mean that any difference between the Monte Carlo programs 
would be less significant. 
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With both two dimensional unfolding and with one dimensional unfolding using 
x COT as the unfolding variable, the spread between the results with the different models 
is reduced compared to one dimensional unfolding results using x vis , and the agree- 
ment between F2GEN and the other models is especially improved. Using different 
unfolding methods with the HERWIG 5.9+/c t (dyn) and PHOJET 1.05 samples makes 
less difference than with the other two Monte Carlo samples, which is as expected for 
models which give a better description of the data. 



As a final test of the unfolding method, several samples of Monte Carlo events 
generated using HERWIG 5.9 with the SaSID structure function were unfolded using 
HERWIG 5.9 with the GRV LO structure function. The results are shown in Figure [17. 
The measured structure function agrees with the input structure function within the 
statistical precision of the Monte Carlo events generated using the SaSID structure 
function. 



The effect of taking into account the shape of f(x) within each bin of x, as in 
Equation |7J, compared to simply extracting as the average input structure function 
of the Monte Carlo sample times the ratio of the number of unfolded events and the 
number of Monte Carlo events in each bin, was also investigated with these samples. 
The effect is largest in the lowest x bins for LEP2 SW, where the decrease in the 
value of K7 when taking into account the shape within the bin is a little less than the 



statistical errors shown in Figure 17 







x 2 






ID 


2D, E^ nt /E tot 


2D, Ef OT /E tot 


(Q 2 ) [GeV 2 ] 


sample 


•2-vis 


•^cor 


•^rec 


•^cor 


•^rec 


•^cor 


1.9 


LEP1 SW 


66.8 


29.1 


22.1 


9.6 


24.7 


26.9 


3.7 




53.1 


26.1 


13.5 


8.1 


19.5 


16.8 


8.9 


LEP1 FD 


15.4 


6.8 


7.1 


3.4 


10.7 


8.7 


17.5 




8.5 


4.5 


8.8 


12.4 


4 


4.1 


10.7 


LEP2 SW 


62.5 


22.4 


18.7 


6.8 


20.6 


16.8 


17.8 




57.2 


18.3 


13 


8.4 


57.7 


15.4 



Table 4: x 2 , as defined in Equation 12, for different unfolding methods. The number 
of bins in x was 4 for the LEP1 SW and LEP2 SW samples, and 3 for the LEP1 FD 
samples. 



6 Results 



The photon structure function is measured by unfolding each data sample in bins 
of log(a;). Each OPAL data sample is divided into two ranges of Q 2 containing ap- 
proximately equal numbers of events. The ranges correspond to (Q 2 ) values of 1.9 and 
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3.7 GeV 2 for the LEP1 SW sample, and 8.9 (10.7) and 17.5 (17.8) GeV 2 for the LEP1 
FD (LEP2 SW) sample, where the average Q 2 values have been obtained from the 
data. F2 is unfolded to a given (Q 2 ) value and does not correspond to an average in 
each bin of Q 2 . 

The results are listed in Table ^. The quoted values were measured as the average 
F2 /a in each bin of x weighted by the unfolded x distribution, according to Equation [7|, 
then corrected to the log centre of each bin, except for the highest x bins where the log 
centre of that portion of the bin below the charm threshold for m c — 1.5 GeV was used. 
The bin-centre corrections are the average of the GRV LO and SaSID predictions for 
the correction from the average F% over the bin to the value of F% at the nominal x 
position. 

The results are also corrected for radiative effects. The radiative corrections were 
calculated using the RADEG [5!J program and are listed in Table [7|, along with the 



bin-centre corrections. The statistical correlations between bins are shown in Table |8|. 

The central value of F^ in each x bin is the average of the data unfolded with 
HERWIG 5.9+/c t (dyn) and PHOJET, using two dimensional unfolding with x cor as the 
first variable and E™ 1 / 'E tot as the second variable. Standard HERWIG 5.9 and F2GEN 
were not used as they are not in acceptable agreement with the data. The systematic 
errors are evaluated by repeating the unfolding with one parameter varied at a time 
and finding the shift in the result. The systematic errors are combined by adding all 
of the individual contributions in quadrature. The summing is done separately for 
positive and negative errors. The systematic effects considered are listed below. 



• Monte Carlo modelling: 

The quoted result is the average of the results obtained using HERWIG 5.9+A; t (dyn) 
and PHOJET 1.05. The errors are symmetrical and equal to half the difference 
between the results with the two Monte Carlo programs. 

• Unfolding method: 

The unfolding was repeated with E ioT /E tot as the second variable, instead of 

/-^tot- 

• Unfolding parameters: 

The number of bins used for the measured variable can be different from the 
number of bins used for the true variable (though it should be at least as large). 
The standard result has 6 bins in the measured variable. This was increased to 
8 to estimate the systematic effects of unfolding. 

• W reconstruction: 

The weighting applied to the energy in the forward calorimeters was varied be- 
tween 2.0 and 3.0. This check allows for uncertainty in the treatment of forward 
energy. 
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• Variations of cuts: 

The composition of the selected events was varied by changing the cuts one at a 
time. The size of the variations reflect the resolution of the measured variables 
and are sufficiently small not to change the average Q 2 of the sample significantly, 
or remove too many events from any single x bin. The variations of the cuts are 
listed in Table ||. 

• Off-momentum electrons: 

There is a possible contamination of off-momentum background in one region of 
azimuthal angle in the 183 GeV data, so as a precaution this <ft region was 
removed in the measurement of K7, and the difference from the result with the 
4> region left in was included as a systematic error. 

• Calibration of the tagging detectors (FD and SW): 

The energy of the tagged electron in the Monte Carlo samples was scaled by 
±1%, to allow for uncertainty in the simulation of the detectors. The size of the 
variation was motivated by a comparison of the E tag / E h distributions in data 
and Monte Carlo. 



• Measurement of the hadronic energy: 

The main uncertainty is in the calibration of the electromagnetic calorimeters 
(excluding the forward region which is dealt with separately). This was varied 
by ±3% in the Monte Carlo samples. The track quality criteria were also varied, 
but the effects were negligible. 



Simulation of background: 

The hadronic background events at LEP1 stem from the production of photons 
or light mesons with a large fraction of the beam energy. The mesons fake 
an electron tag due to their decays into photons. The cross-section for these 
events was measured by OPAL with an accuracy of about 50%, and found to be 
consistent with the JETSET prediction p4 |. The normalisation of the simulated 



background was therefore varied by ±50%. 



• Bin-centre correction: 

The corrections depends on the shape of in each bin. The average of the 
corrections based on GRV LO and SaS ID was used, as these parameterisations 
are the closest to the data. The error is half the difference between the GRV LO 
and SaS ID corrections, and is symmetric. 



Because the estimation of each source of systematic error has some statistical fluc- 
tuation due to the changing distributions (e.g. when a cut is varied), the quadratic sum 
of all individual sources is an overestimate of the total systematic error. The expected 
statistical component of each source of systematic error that changes the data sample 
is determined by unfolding 8 Monte Carlo samples, each about the same size as the 
data sample. 
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The corrected systematic error on the ith bin from source k is then 



(A'f hk ) 2 = m,kf - cr 2 Akk for Af 2 k > a\ h k 

^'fi,k = otherwise (13) 

where Afc^ is the shift on the ith bin in the data when changing a parameter k, and 
o"A/i k is the expected statistical component of the shift, which is approximated by the 
statistical spread in the systematic error estimates for source k in the 8 Monte Carlo 
samples. In a few cases this was larger than the statistical error in the data. In these 
cases o"A/ lfc was set to the statistical error in the data, in order not to hide possibly 
significant systematic errors. This procedure means that for the individual systematic 
errors either the expected statistical component is subtracted, or the error is set to zero 
if the observed shift in the data is consistent with a statistical fluctuation as predicted 
by the 8 Monte Carlo samples. 

The total systematic error is the quadratic sum of all the individual contributions 
A'fi t k from the above sources. The sum is made separately for deviations above and 
below the standard result. The systematic errors from each source as a function of x 
and Q 2 are listed in Table 9. 



The results are shown along with previous OPAL measurements of i 7 ^ [16] in Fig- 
ures [18| and [1^. The overlapping results from the LEP1 FD and LEP2 SW samples are 
in good agreement. The measurements of F% using LEP1 data for (Q 2 ) = 1.9 GeV 2 
and (Q 2 ) = 3.8 GeV 2 are lower than the previous OPAL LEP1 results, which were un- 
folded using HERWIG 5.8d. Repeating the unfolding with HERWIG 5.8d gives results 
which are consistent with the old analysis, but with better precision. The HERWIG 
5.8d Monte Carlo model has now been replaced by HERWIG 5.9+/c t (dyn), which gives 
a better description of the data. The results in the four higher Q 2 bins are consis- 
tent with previous OPAL measurements, and have smaller errors. The previous LEP1 
results with electrons tagged in SW or FD are superseded by the present analysis. 

In Figures |2(] and ^1] the results are compared to measurements of F% from other 
experiments: TPC/2 7 0, PLUTO 0, TOPAZ p|, ALEPH p|, DELPHI |7| and 



L3 0. Also shown in Figures gU| and [H| are the GRV LO, SaSID and WHIT1 §|] 
parameterisations of F% , and the naive quark-parton model (QPM). The QPM predic- 
tion, which only models the point-like component of i 7 ^, is calculated for four active 
flavours with masses of 0.2 GeV for light quarks and 1.5 GeV for charm quarks. 

The previous results are found to be generally consistent with the new OPAL mea- 
surements. The largest differences are from the older measurement by TPC/27 at 
(Q 2 ) = 2.8 GeV 2 , which suggests a different shape of to all the other measure- 
ments. The L3 results obtained at (Q 2 ) = 5.0 GeV 2 are consistently higher than the 
OPAL measurements. However, because the L3 points are highly correlated, due to 
the finer binning in x, the discrepancy looks stronger than it actually is. 

7 The starting scale of the evolution for the WHIT parameterisation is 4 GeV 2 ; consequently it can 
only be compared to the measurements above (Q 2 ) — 4 GeV 2 . 
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cut 


LEP1 SW 


LEP1 FD 


LEP2 SW 


E tag /E h min 


±0.025 (0.75) 


±0.025 (0.775) 


#tag mm [mrad] 


+2 (27) 


+2 (60) 


+2 (33.25) 


^tag max [mrad] 


-2 (55) 


-2 (120) 


-2 (55) 


EjE h max 


±0.05 (0.25) 


W vis min [GeV] 


+1.0 (2.5) 


W vis max [GeV] 


±5 (40) 


±5 (60) 



Table 5: Systematic variations in the cuts. The standard cuts are given in brackets 
after the variations. 

Due to large spread in the theoretical predictions for the P 2 suppression of F% (as 
discussed in Section f|) the OPAL measurements are not corrected for this effect. Since 
the data contain a P 2 suppression, which is not included in the Monte Carlo simulation, 
applying the correction would lead to an x dependent increase of the measured K7. This 
should be taken into account when comparing to the parameterisations of K7, which 
are all shown for P 2 = 0. In general the shape of the GRV LO parameterisation is 
consistent with the OPAL data in all the accessible x and Q 2 regions. The normalisation 
is also consistent with the data, except at the lowest scale, (Q 2 ) = 1.9 GeV 2 , where 
GRV is too low. The SaSID LO prediction shows a slower evolution with Q 2 than the 
GRV prediction. At low Q 2 it is slightly above GRV LO, whereas at the largest Q 2 
values shown, it falls below GRV LO. Within the precision of the OPAL measurement, 
the description of the data by SaSID LO is of similar quality to GRV LO. In the region 
of applicability, the WHIT1 prediction is higher than the OPAL data and flatter than 
the other predictions, though the shape is still consistent with the data. 

The hadron-like component is predicted to dominate at low values of x. For x < 0.1 
the naive quark-parton model is not able to describe the data, indicating that the 
photon must contain a significant hadron-like component at low x. 



7 Conclusions 

The photon structure function F% has been measured using deep inelastic electron- 
photon scattering events recorded by the OPAL detector during the years 1993-1995, 
1997 and 1998, at e+e" centre-of-mass energies of 91 GeV (LEP1), and 183-189 GeV 
(LEP2). F2 has been measured as a function of x to the lowest attainable x values, in 
six ranges of Q 2 (including two overlapping pairs) corresponding to average Q 2 values 
of 1.9, 3.7 GeV 2 for LEP1 SW, and 8.9 (10.7), 17.5 (17.8) GeV 2 for LEP1 FD (LEP2 
SW). In previous OPAL studies of the photon structure function, it became clear that 
a large source of uncertainty in the measurement came from the Monte Carlo modelling 
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of the hadronic final state of deep inelastic electron-photon scattering events. Since 
then, improved Monte Carlo models have become available. In a comparison of the 
energy flows and E^ ut distributions, which are very sensitive to the modelling of the 
hadronic final state, these new models, HERWIG 5.9+£; t (dyn) and PHOJET 1.05, 
give better descriptions of OPAL data than the HERWIG 5.9 and F2GEN programs. 
Consequently the latter two programs have not been used for the K7 measurement. 
Previous OPAL measurements of using LEP1 data with electrons tagged in SW or 
FD are superseded by this analysis. 

To further reduce the Monte Carlo modelling error, two dimensional unfolding has 
been introduced, using / E tot as a second unfolding variable. Also, the reconstruc- 
tion of the invariant mass of the hadronic final state has been improved by including 
information from the deeply inelastically scattered electron, and by scaling the energy 
observed in the forward calorimeters to partially compensate for energy losses. Monte 
Carlo modelling of the final state is still a significant source of systematic error, but it 
no longer dominates all other sources. The total systematic errors are of comparable 
size to the statistical errors. 

Although the precision of the measurement at low x has been considerably improved 
it is still insufficient to determine whether or not there is a rise in F^ in that region. 

The GRV LO and SaSID parameterisations are generally consistent with the OPAL 
data in all the accessible x and Q 2 regions, with the exception of the measurement at 
the lowest scale, (Q 2 ) = 1.9 GeV 2 , where GRV is too low. In contrast, the naive quark- 
parton model is not able to describe the data for x < 0.1. These results show that the 
photon must contain a significant hadron-like component at low x. 
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(Q 2 > [GeV 2 ] 


sample 


bin 


x range 


X 


F 2 > 


1.9 


LEP1 SW 


I 

II 
III 

IV 


0.0006 < x < 0.0028 
0.0028 < x < 0.0143 
0.0143 < x < 0.0724 
0.0724 < x < 0.3679 


0.0012 
0.0063 
0.0322 
0.1124 


0.269 ±0.0271^34 
0.177 ±0.009lHiI 
0.179 ± 0.007 t° °o l 
0.227 =b 0.010±g;gil 


3.7 


LEP1 SW 


I 

II 
III 

IV 


0.0015 < x < 0.0067 
0.0067 < x < 0.0302 
0.0302 < x < 0.1353 
0.1353 < x < 0.6065 


0.0032 
0.0143 
0.0639 
0.1986 


0.269 ±0.0331^33 
0.232 ±0.0131^21 
0.259 ± 0.010±g;8Si 
0.296 ±0.0141^1 


8.9 


LEP1 FD 


I 

II 
III 


0.0111 < x < 0.0498 
0.0498 < x < 0.2231 
0.2231 < x < 0.8187 


0.0235 
0.1054 
0.3331 


0.221 ±0.0171^26 
0.308 ± 0.014 ±%%\l 
0.379 ±0.0221^15 


10.7 


LEP2 SW 


I 

II 
III 

IV 


0.0009 < x < 0.0050 
0.0050 < x < 0.0273 
0.0273 < x < 0.1496 
0.1496 < x < 0.8187 


0.0021 
0.0117 
0.0639 
0.3143 


0.362 ± 0.045 ±g;g|| 
0.263 ± 0.015 ±8:Slo 
0.275 ± 0.011 ±g;g|g 

o.35i ± 0.012 ±8:gfi 


17.5 


LEP1 FD 


I 

II 
III 


0.0235 < x < 0.0821 
0.0821 < x < 0.2865 
0.2865 < x < 0.9048 


0.0439 
0.1534 
0.3945 


0.273 ±0.028l8;839 
0.375 ± 0.023 18:o?3 

o.5oi ± 0.027 18:8?9 


17.8 


LEP2 SW 


I 

II 
III 

IV 


0.0015 < x < 0.0074 
0.0074 < x < 0.0369 
0.0369 < x < 0.1827 
0.1827 < x < 0.9048 


0.0033 
0.0166 
0.0821 
0.3483 


0.428 ± 0.061 18:o7i 
0.295 ± 0.019 l8:oi 
0.336 ±0.01318:842 
0.430 ± 0.013 18:o25 



Table 6: Results for K7 /a as a function of x for four active flavours in bins of Q 2 . 
The first errors are statistical and the second systematic. The structure function was 
unfolded in bins defined by the x ranges and corrected to the x values given. 
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(Q 2 ) [GeV 2 ] 


sample 


bin 


x range 


X 


radiative 
correction 


bin-centre 
correction 


1.9 


LEP1 SW 


I 


0.0006 < x < 0.0028 


0.0012 


-12.7 


-4.2 






II 


0.0028 < x < 0.0143 


0.0063 


-9.0 


0.4 






III 


0.0143 < x < 0.0724 


0.0321 


-7.1 


1.8 






IV 


0.0724 < x < 0.3679 


0.1124 


-6.0 


4.7 


3.7 


LEP1 SW 


I 


0.0015 < x < 0.0067 


0.0032 


-11.8 


-5.0 






II 


0.0067 < x < 0.0302 


0.0143 


-8.9 


0.6 






III 


0.0302 < x < 0.1353 


0.0639 


-7.3 


1.9 






IV 


0.1353 < x < 0.6065 


0.1986 


-6.5 


1.6 


8.9 


LEP1 FD 


I 


0.0111 < x < 0.0498 


0.0235 


-7.7 


0.9 






II 


0.0498 < x < 0.2231 


0.1054 


-6.3 


2.4 






III 


0.2231 < x < 0.8187 


0.3331 


-4.1 


-0.9 


10.7 


LEP2 SW 


I 


0.0009 < x < 0.0050 


0.0021 


-12.5 


-8.7 






II 


0.0050 < x < 0.0273 


0.0117 


-7.3 


-0.5 






III 


0.0273 < x < 0.1496 


0.0639 


-4.4 


3.2 






IV 


0.1496 < x < 0.8187 


0.3143 


-2.2 


-1.0 


17.5 


LEP1 FD 


I 


0.0235 < x < 0.0821 


0.0439 


-9.4 


2.1 






II 


0.0821 < x < 0.2865 


0.1534 


-7.9 


2.5 






III 


0.2865 < x < 0.9048 


0.3945 


-6.5 


0.0 


17.8 


LEP2 SW 


I 


0.0015 < x < 0.0074 


0.0033 


-13.6 


-8.2 






II 


0.0074 < x < 0.0369 


0.0166 


-9.9 


-0.5 






III 


0.0369 < x < 0.1827 


0.0821 


-8.4 


3.7 






IV 


0.1827 < x < 0.9048 


0.3483 


-7.3 


-0.4 



Table 7: Corrections to the result as a function of x in bins of Q 2 , as a percentage 
of the non-corrected F^. The radiative corrections were predicted by RADEG @. 
The bin-centre corrections are the average of the GRV LO and SaSID predictions for 
the correction from the average over the bin to the value of at the nominal x 
position. The x positions are at the log centre of the bins, except for the highest x 
bins, where they are at the the log centre of that portion of the bin below the charm 
threshold for m c = 1.5 GeV. 
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(Q 2 )=1.9 GeV 2 (LEP1 SW) 





I 


II 


III 


IV 


I 


1.00 








II 


-0.28 


1.00 






III 


0.03 


-0.35 


1.00 




IV 


0.01 


0.10 


-0.48 


1.00 



(Q 2 )=8.9 GeV 2 (LEP1 FD) 





I 


II 


III 


I 


1.00 






II 


-0.02 


1.00 




III 


-0.12 


-0.31 


1.00 



(Q 2 ) = 17.5 GeV 2 (LEP1 FD) 





I 


II 


III 


I 


1.00 






II 


-0.09 


1.00 




III 


-0.09 


-0.35 


1.00 



Table 8: Statistical correlations between 
the bins listed in Table |6|. 



(Q 2 )=3.7 GeV 2 (LEP1 SW) 





I 


II 


III 


IV 


I 


1.00 








II 


-0.27 


1.00 






III 


0.02 


-0.34 


1.00 




IV 


0.01 


0.10 


-0.52 


1.00 



(Q 2 )=10.7 GeV 2 (LEP2 SW) 





I 


II 


III 


IV 


I 


1.00 








II 


-0.32 


1.00 






III 


0.04 


-0.29 


1.00 




IV 


0.00 


0.06 


-0.36 


1.00 



(Q 2 )=17.8 GeV 2 (LEP2 SW) 





I 


II 


III 


IV 


I 


1.00 








II 


-0.35 


1.00 






III 


0.05 


-0.33 


1.00 




IV 


-0.01 


0.09 


-0.41 


1.00 



for each sample. The numerals refer to 
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3.62 
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0.02 


0.45 
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LEP1 
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17 3 
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-4 23 








-1 62 


o 










-9 83 


-1 42 


4 7 








4 06 




4 53 


2 77 


-3 55 


1 36 








0143 


232 






9 2 


-8 29 


4 28 


o 


o 








o 


-2 43 


98 




22 


-1 47 






-2 18 


2 69 


-1 41 


1 22 


12 


-0 27 


35 








0.0639 


0.259 


3.9 


2.4 


4.9 


-0.5 


-0.84 


-3.35 


o 


o 


o 


-0.92 


o 


0.15 


o 


-0.36 


-0.07 


0.6 


-0.16 




-2.35 


1.94 


-2.37 


1.25 


-0.01 


0.04 


0.1 








0. 1986 


0.296 


4.6 


9.9 


7.3 


6.59 


-1.68 


6.01 


-1.01 


o 


o 


0.32 


Q 


o 


o 


0.17 


0.01 


-0.26 


o 




-1.67 


2.21 


3.71 


-1.71 


-0.01 


-0.02 


0.54 


8.9 


LEP1 


FD 


0.0235 


0.221 


7.6 


13.5 


11.6 


-10.51 


o 


3.85 


6.41 


1.95 


-1.65 


o 


o 


o 


-0.96 


o 


0.01 


-1.02 


o 


— 


-2.12 


1.29 


-2.68 


2.31 


2.14 


-2.76 


0.84 








0.1054 


0.308 


4.6 


3.6 


3.9 


0.64 


-2.68 








-2.05 








2.47 





0.53 
















-0.42 


0.62 


-1.75 


2.29 


-0.36 


0.46 


0.06 








0.3331 


0.379 


5.8 


4.5 


4 


-2.63 


1.67 








1.04 


























0.08 




-0.24 


-0.74 


3.04 


-2.91 


0.13 


-0.22 


0.37 


10.7 


LEP2 


SW 


0.0021 


0.362 


12.3 


16.1 


10.7 


-1.21 








4.23 


-3.94 


-5.05 


9.62 


3.5 











3.45 





4.08 


7.41 


-8.11 


6.85 


-2.34 


0.27 








1.54 








0.0117 


0.263 


5.8 


12.3 


11.6 


-9.85 


5.28 


3.68 














-5.96 





0.35 





-0.4 








3.27 


0.67 


-0.22 


-0.86 


1.3 








0.5 








0.0639 


0.275 
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10.4 


10.8 


-9.14 





-4.5 








-2.17 


0.64 


4.44 





-2.32 





0.09 








1.14 


0.83 


-0.45 


-1.5 


1.64 








0.41 








0.3143 


0.351 


3.5 
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4.5 


-3.68 




















2.74 





3.87 





-0.05 


0.02 





2.57 


1.74 


-1.83 


1.76 


-1.87 








0.17 


17.5 


LEP1 


FD 


0.0439 


0.273 


10.4 


11.8 


14.4 


-8.93 











-1.21 


-3.38 








2.55 


0.95 


-0.32 


-0.03 





-6.56 





-3.43 


2.75 
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3.87 


5.23 


-7.45 


0.96 








0.1534 


0.375 


6.1 


5.4 


3.6 


-1.34 





3.25 








0.52 


1.26 





0.66 





0.23 


-0.02 


-0.92 


1.25 
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0.03 
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-0.02 











-2.72 


2.76 
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-1.61 


-0.05 
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1.42 


17.8 


LEP2 


SW 


0.0033 


0.428 


14.2 


12.8 


16.7 


-2.22 


-2.88 


-2.91 








-15.09 


2.33 








-0.57 








2.89 





0.68 


-5.28 


11.84 


-1.03 


1.26 








1.33 








0.0166 


0.295 


6.4 


11.2 


6.8 


-5.79 


6.79 





2.84 





2.04 


1.95 





-0.78 


1.03 





0.52 








4.49 
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2.63 


-0.99 


0.99 








0.81 








0.0821 
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3.9 


12.3 


12.6 


-11.56 


-1.87 





-0.19 





-1.68 
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-0.48 


0.43 


-1.05 


-2.41 


3.32 
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0.18 





-3.29 


5.65 





0.02 
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-2.47 
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Table 9: Systematic errors for F^jot as a function of x in bins of Q 2 . The statistical error and the contributions to the systematic 
error from each of the sources are listed as percentages of F% /a . The systematic errors are corrected for the expected statistical 
component as described in the text. They are summed in quadrature separately for positive and negative deviations, except for 
'MC modelling' and 'bin centre' which are symmetric. 



Vl 


1600 






> 


1400 




1200 




1000 




800 




600 




400 




200 








~i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 

• OPALLEP1SW (a) 

— HERWIG 5.9+k t (dyn) 

— - PHOJET 1.05 

" H y*y->x + x~ 

— total background 




0.7 0.8 0.9 




VI 






1400 


> 








12UU 




1 AAA 




QAA 




600 




400 




zuu 







VI 




a 4500 






£ 4000 




3500 




3000 




2500 




2000 




1500 




1000 




500 









40 50 

tag t mrad ] 




Q 2 [GeV 2 ] 



Figure 2: Comparison of data distributions with Monte Carlo predictions for the LEP1 
SW sample. The dominant background source, 7*7 — > t + t~ , the total background 
and the sum of the signal and the total background for HERWIG 5.9+fc t (dyn) and 
PHOJET 1.05 are shown. The Monte Carlo samples have been normalised to the data 
luminosity. All selection cuts have been applied, except for any cut on the variable in 
the plot. The cuts are shown as dotted lines. The distributions shown are: a) E tllg /E h , 
the energy of the tagged electron as a fraction of the beam energy, b) tag , the polar 
angle of the tagged electron, c) the measured Q 2 and d) E a /Eb, the energy of the most 
energetic electromagnetic cluster in the hemisphere opposite the tagged electron, as a 
fraction of the beam energy. The errors are statistical only. 
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Figure 3: Comparison of data distributions with Monte Carlo predictions for the LEP1 
FD sample. The dominant background sources, 7*7 — > t + t~ and Z° — > hadrons, the 
total background and the sum of the signal and the total background for HERWIG 
5.9+/c t (dyn) and PHOJET 1.05 are shown. The Monte Carlo samples have been nor- 
malised to the data luminosity. All selection cuts have been applied, except for any 
cut on the variable in the plot. The cuts are shown as dotted lines. The variables in 
the plots are as defined in Figure |[ The errors are statistical only. 
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Figure 4: Comparison of data distributions with Monte Carlo predictions for the LEP2 
SW sample. The dominant background source, 7*7 — > r + r~, the total background 
and the sum of the signal and the total background for HERWIG 5.9+/c t (dyn) and 
PHOJET 1.05 are shown. The Monte Carlo samples have been normalised to the data 
luminosity. All selection cuts have been applied, except for any cut on the variable in 
the plot. The cuts are shown as dotted lines. The variables in the plots are as defined 
in Figure [|. The errors are statistical only. 
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Figure 5: Comparison of data distributions with Monte Carlo predictions for the 
LEP1 SW sample. The dominant background sources, 7*7 — > t + t~, Z° — > hadrons 
and 7*7* — > hadrons (with P 2 >1.0 GeV 2 ), the total background and the sum of the 
signal and the total background for HERWIG 5.9+A; t (dyn) and PHOJET 1.05 are 
shown. For illustration, the background from 7*7 — * yU + /i~ is also shown, though it is 
only significant below the iV ch cut. The Monte Carlo samples have been normalised 
to the data luminosity. All selection cuts have been applied, except for any cut on 
the variable in the plot. The cuts are shown as dotted lines. The distributions shown 
are: a) W v i S , the measured invariant mass of the hadronic final state, b) log(x v is), 
the logarithm of the measured value of x. c) iV c h, the number of tracks in the event 
and d) E^, the transverse component of energy out of the tag plane. The errors are 
statistical only. 
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Figure 6: Comparison of data distributions with Monte Carlo predictions for the 
LEP1 FD sample. The dominant background sources, 7*7 — > t + t~, Z° — > hadrons 
and 7*7* — >• hadrons (with P 2 >1.0 GeV 2 ), the total background and the sum of the 
signal and the total background for HERWIG 5.9+£; t (dyn) and PHOJET 1.05 are 
shown. For illustration, the background from 7*7 — > / u + fi" is also shown, though it is 
only significant below the cut. The Monte Carlo samples have been normalised to 
the data luminosity. All selection cuts have been applied, except for any cut on the 
variable in the plot. The cuts are shown as dotted lines. The variables in the plots are 
as defined in Figure |5|. The errors are statistical only. 
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Figure 7: Comparison of data distributions with Monte Carlo predictions for the 
LEP2 SW sample. The dominant background sources, 7*7 — > t + t~ , Z° — > hadrons 
and 7*7* — > hadrons (with P 2 >4.5 GeV 2 ), the total background and the sum of the 
signal and the total background for HERWIG 5.9+£; t (dyn) and PHOJET 1.05 are 
shown. For illustration, the background from 7*7 — > / u + fi" is also shown, though it is 
only significant below the cut. The Monte Carlo samples have been normalised to 
the data luminosity. All selection cuts have been applied, except for any cut on the 
variable in the plot. The cuts are shown as dotted lines. The variables in the plots are 
as defined in Figure |5|. The errors are statistical only. 
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Figure 8: Comparison of hadronic energy flow per event for data and Monte Carlo 
as a function of the pseudorapidity, r\ = — log(tan(#/2)), where the polar angle 6 is 
measured with respect to the beam axis on the side of the tagged electron. All three 
samples, LEP1 SW, LEP1 FD and LEP2 SW are shown. The tagged electron is not 
included in these plots. The errors are statistical only. 
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Figure 9: Data and signal Monte Carlo distributions for two alternative variables for 
two dimensional unfolding, for the LEP1 SW sample divided into three bins of x v i S . 
Plots a) c) and e) show E{ OT /E tot , the observed energy in the forward region divided by 
the total observed energy. Plots b) d) and f) show E™ 1 / E tot , defined as the transverse 
component of hadronic energy out of the plane of the tagged electron, divided by the 
total energy. The Monte Carlo samples have been normalised to the number of data 
events in each plot. The errors are statistical only. 
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Figure 10: Data and signal Monte Carlo distributions for two alternative variables for 
two dimensional unfolding, for the LEP1 FD sample divided into three bins of x v ; s . The 
variables are as defined in Figure The Monte Carlo samples have been normalised 
to the number of data events in each plot. The errors are statistical only. 
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Figure 11: Data and signal Monte Carlo distributions for two alternative variables for 
two dimensional unfolding, for the LEP2 SW sample divided into three bins of x v i S . The 
variables are as defined in Figure The Monte Carlo samples have been normalised 
to the number of data events in each plot. The errors are statistical only. 
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Figure 12: Tests of the unfolding procedure using the LEP1 SW low Q 2 data unfolded 
with the HERWIG 5.9+A; t (dyn) Monte Carlo sample, (a) Comparison of the results 
using GURU, RUN and BAYES. (b) Different degrees of freedom for one dimensional 
unfolding, (c) Different measured x variables for one dimensional unfolding, (d) Dif- 
ferent degrees of freedom for two dimensional unfolding using a random number as the 
second variable. The errors are statistical only. 
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Figure 13: (a) The correlation between the generated and measured invariant mass, W 
and ^measured, for the HERWIG 5.9+£; t (dyn) LEP2 SW Monte Carlo sample, using, 
Wvis, W rec and W cor . (b) The corrected energy, E COT) (2.5 times the visible energy) in 
the forward region against the total generated energy, E, deposited in that region, for 
HERWIG 5.9+/ct(dyn) and PHOJET 1.05. The vertical error bars represent the RMS 
spread within each bin. The dotted lines represent perfect correlation. 
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Figure 14: Unfolding of the LEP1 SW low Q 2 sample with 4 Monte Carlo models 
and different unfolding variables, a) x v i S , b) x cor , c) x rcc and E^/Etot, d) rc cor and 
E^ ut /E tot , e) x re c and E ior /E tot , f) rr cor and E ior /E tot . The errors are statistical only. 
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Figure 15: Unfolding of the LEP1 FD low Q 2 sample with 4 Monte Carlo models and 
different unfolding variables. The sub-figures are as described in Figure [14]. The errors 
are statistical only. 
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Figure 16: Unfolding of the LEP2 SW low Q 2 sample with 4 Monte Carlo models and 
different unfolding variables. The sub-figures are as described in Figure [14]. The errors 
are statistical only. 
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Figure 17: Two independent samples of Monte Carlo events generated using HERWIG 
5.9 with the SaSID structure function unfolded with HERWIG 5.9 with the GRV LO 
structure function, in six Q 2 ranges. The solid histogram shows the SaS ID structure 
function at the same average Q 2 as the sample in each plot, weighted by the x distri- 
bution in the HERWIG SaSID sample (this is the quantity measured by Equation |7|). 
The dotted lines show the GRV LO structure function weighted by the unfolded x dis- 
tribution. The unfolded SaSID samples were about the same size as the data samples 
in each Q 2 region. The errors are statistical only. 
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Figure 18: The measurement of F^/a using the LEP1 SW sample, for (Q 2 ) values 
of (a) 1.9 and (b) 3.7 GeV 2 . Also shown are the previous OPAL results in these 
Q 2 ranges, which were unfolded using HERWIG 5.8d, and the result of unfolding the 
LEP1 SW data using HERWIG 5.8d. For each point, the inner error bars show the 
statistical error and the full error bars show the total error, except for the new result 
with HERWIG 5.8d, for which only statistical errors are shown. The positions of the 
new OPAL points are as given in Table [| The other points are shown at the centre of 
the log(x) bin. The curves show the GRV LO structure function. 
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Figure 19: The measurement of ja using the LEP1 FD and LEP2 SW samples for 
(Q 2 ) values of (a) 8.9 (10.7) and (b) 17.5 (17.8) GeV 2 for LEP1 (LEP2). Also shown 
are the previous OPAL results in these Q 2 ranges, which were unfolded with HERWIG 
5.8d ((Q 2 >=7.5 GeV 2 and (Q 2 )=14.7 GeV 2 ) and HERWIG 5.9 ((Q 2 )=9.0 GeV 2 and 
(Q 2 )=14.5 GeV 2 ) using a linear x scale. For each point, the inner error bars show the 
statistical error and the full error bars show the total error. The positions of the new 
OPAL points are as given in Table |6|. The other points with closed symbols are shown 
at the centre of the log(x) bin, and those with open symbols are shown at the average 
x value of the bin. The curves show the GRV LO structure function. 
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Figure 20: The measurement of F^/a using the LEP1 SW sample, for (Q 2 ) values of 
(a) 1.9 and (b) 3.7 GeV 2 . Also shown are the results from L3 |18] PLUTO [10J, and 
TPC/27 For L3 the two sets of points were unfolded using different Monte Carlo 
programs. The lower / upper points correspond to PHOJET 1.05 / TWOGAM. For 
each point, the inner error bars show the statistical error and the full error bars show 
the total error. The positions of the new OPAL points are as given in Table |6|. The 
other points with closed symbols are shown at the centre of the log (a;) bin, and those 
with open symbols are shown at the average x value of the bin. The curves show the 
GRV LO, SaSID, WHIT1 and QPM structure functions. 
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Figure 21: The measurement of F^/oc using the LEP1 FD and LEP2 SW samples 
for (Q 2 ) values of (a) 8.9 (10.7) and (b) 17.5 (17.8) GeV 2 for LEP1 (LEP2). Also 
shown is a selection of results from other experiments: ALEPH ||19|| , DELPHI (17 



L3 ||18||, PLUTO [TO], and TOPAZ ||15||. For each point, the inner error bars show the 



statistical error and the full error bars show the total error. The positions of the new 
OPAL points are as given in Table |6|. The other points with closed symbols are shown 
at the centre of the log(x) bin, and those with open symbols are shown at the average 
x value of the bin. The curves show the GRV LO, SaSID, WHIT1 and QPM structure 
functions. 
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